Interlayer Coupling and p-wave Pairing in Strontium Ruthenate 
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On the basis of a three orbital model and an effective attractive interaction between electrons we 
investigate the possible superconducting states, with p and /-wave internal symmetry, of S^RuCU. 
For an orbital dependent interaction which acts between in plane and out of plane nearest neighbour 
Ruthenium atoms we find a state for which the gap in the quasi-particle spectra has a line node on 
the a and j3 sheets of the Fermi Surface, but it is complex with no nodes on the 7-sheet. We show 
that this state is consistent with all the available experimental data. In particular, we present the 
results of our calculations of the specific heat and penetration depth as functions of the temperature. 

PACS numbers: 74.70.Pq, 74.20.Rp, 74.25.Bt 



Ever since triplet pairing was discovered in superfluid 
3 He Q, during the early seventies, there has been a con- 
stant search for the superconducting analogue of this in- 
triguing macroscopic quantum phenomenon. Although, 
as yet, there is no metal for which triplet paring has been 
definitively demonstrated there are a number of good 
candidates. The evidence that Sr2Ru04 is a triplet su- 
perconductor is particularly strong ||||. Nevertheless, 
even in this case the full symmetry of the equilibrium 
state below T c remains open to debate P~^2[. 

One of the puzzles currently at the centre of atten- 
tion is the apparent incompatibility between experimen- 
tal evidence for broken time-reversal symmetry in the su- 
perconducting state [T^JTo[ | and equally convincing mea- 
surements indicating that the order parameter d(k) has 
a line of nodes on the Fermi Surface pT|Jl5| ]. The rea- 
son why this state of affairs represents a dilemma is that 
for all odd parity spin triplet pairing states in tetrago- 
nal crystals, group theory does not require the simulta- 
neous presence of both broken time-reversal symmetry 
and line nodes |l6[ ). Consequently, due to their lower 
condensation energy, line nodes are unlikely. For in- 
stance, the pairing state d(k) ~ (k x + ik y )e z , proposed 
by Agterberg et al. [Q on the grounds that it minimizes 
the free energy, obviously breaks time reversal invari- 
ance and has no line nodes. Of course, such states as 
d(k) ~ (k x + iky)(k x k y )e z discussed by Graf and Bal- 
atsky |Q and other f-wave states |^|,^ 11 are allowed by 
symmetry considerations but, the point is that, the nodes 
are not required by symmetry and hence it is not very 
attractive ansatz to build a theory on. Under these cir- 
cumstances it is more advantageous to study physically 
motivated microscopic models even if the question of the 
actual mechanism of pairing is to be avoided. In this let- 
ter we propose and investigate one such physically moti- 
vated model. 



Our model is prompted by the observation, of 
Hasegawa et al. that coupling between the Ruthe- 
nium layers in Sr2Ru04 leads to convenient, horizontal, 
k z = ±- lines of zeros on the Fermi Surface. It features 
two, intra and inter plane, interaction constants U\\ and 
U± respectively, which describe the attraction between 
electrons each occupying one of three, ta ff , orbitals on 
Ru atoms which are nearest neighbours either in plane 
and out of plane. As we shall show, this simple physi- 
cal picture yields horizontal nodes on the a, (3 sheets of 
the Fermi Surface whilst the 7-sheet is fully gaped with 
d(k) ~ (sinfca; + isink y )e z , in quantitative as well as 
qualitative agreement with a number of experiments. 

Hasegawa et al. |ll| , treated the case of a single band 
only and, unlike us, they made no quantitative contact 
with experiments. The more recent work of Zhitomirsky 
and Rice jl2| is closer to ours, although dealing with a 
two band model only. They also couple electrons on dif- 
ferent layers and find line nodes on the a, (3 sheets but 
a gap on 7. However, the physics behind their model, 
as well as its consequences, are quite different from ours. 
Largely, this is due to a difference in strategy. In a multi- 
band BCS like model, with different coupling constants 
for each band, one generically finds multiple phase tran- 
sitions as the different sheets of the Fermi Surface are 
gaped on lowering the temperature. Since experimentally 
there is only one jump in the specific heat, at T c = 1.5K, 
in constructing a sensible model one must tailor it to 
eliminate such double transitions. Zhitomirsky and Rice 
[fiJI chose to couple the order parameters on different 
sheets of the Fermi Surface. This hybridized the differ- 
ent order parameters and led to a single transition. They 
referred to such hybridization as a inter-band proximity 
effect. We, on the other hand, adjusted the two cou- 
pling constants, Un and U±, so that the transition on 
the a, (3 sheets would occur at more or less the same 
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temperature as as on 7. Somewhat surprisingly these 
two approaches imply different physics. The proximity 
model requires a three point interaction to mix the in- 
plane and out of plane Cooper pairs whilst our "local 
bond" model is a strictly a two point effect. Moreover, 
the "local bond" model has fewer free parameters by con- 
struction and therefore, as will be shown below, is more 
readily compared with experiments. 

To describe the superconducting state we employ a 
simple multi-band attractive Hubbard model: 

H ^ ^ ((^m f^)^ij &mm' ^mm'(^j)) ^merA? m/(T 

ijmm' ,a 

-\ U°° m ,(ii)n lma h jm > al (1) 

ijmm' a a' 

where m and ml refer to the three Ruthenium t 2g orbitals 
a = xz, b — yz and c — xy and i and j label the sites of a 
body centered tetragonal lattice. The hopping integrals 
tmm'(ij) and site energies e m were fitted to reproduce 
the experimentally determined Fermi Surface |l7| . The 
set of interaction constants U™ m , (ij) describe attraction 
between electrons on nearest neighbour sites with spins a 
and a' and in orbitals m and ml . Thus our actual calcu- 
lations consists of solving, self-consistently, the following 
Bogoliubov-de Gennes equation: 
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where H mm '(ij) is the normal spin independent part of 
the Hamiltonian, and the AJ^L, (ij) is self consistently 
given in terms of the pairing amplitude, or order param- 
eter C(y'), 
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defined by the usual relation 
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where v enumerates the solutions of Eq. Q. 

We solved the above system of Bogoliubov de Gennes 
equations including all three bands and the experimen- 
tal three dimensional Fermi surface. We assumed that 
the pairing interaction U™ m i(ij) for nearest neighbours 
in plane is only acting for the c (d xy ) Ru orbitals. We 
further assumed that the nearest neighbour inter plane 
interaction acts only in a and b orbitals (d xz , d yz ). The 
motivation for this is that the dominant hopping inte- 
grals in plane are between c orbitals, and the largest out 
of plane hopping integrals are for a and b. Therefore we 
have only two coupling constants C/|| and U± describing 
these two physically different interactions. Our strat- 
egy is to adjust these phenomenological parameters in 



order to obtain one transition at the experimentally de- 
termined T c . Thus, beyond fitting T c , there are no fur- 
ther adjustable parameters, and one can compare directly 
the calculated physical properties of the superconducting 
states to those experimentally observed. Consequently, 
if one obtains a good overall agreement one can say that 
one has empirically determined the form of the pairing 
interaction in a physically transparent manner. 
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FIG. 1. Order parameters, \A x aa \, \A%\, |A* a |, |A£„| as 
functions of temperature, (dashed lines), and excluding z and 
/ (full lines). 

Because the pairing interactions U™ m ,(ij) were as- 
sumed to act only for nearest neighbour sites in or 
out of plane, the pairing potential A^ m , (ij) is also re- 
stricted to nearest neighbours. We further focus on 
only odd parity (spin triplet) pairing states for which 
the vector d ~ (0,0, i.e. A^ m ,(ij) = A# TO ,(y), 
and AH m ,(ij) = A^ m ,(ij) = 0. Therefore in gen- 
eral we have the following non-zero order parameters 
(i) for in plane bonds: A cc (e x ), A cc (ej,), and (ii) for 
inter-plane bonds: A ao (R y ), A^Ry), A bb (R 2j ) for 
TLij = (±a/2,±a/2,c/2). 

Taking the lattice Fourier transform of Eq. H] the cor- 
responding pairing potentials in k-space have the general 
form (suppressing the spin indices for clarity): 



A cc (k) = A x cc sin k x + A v cc sin k y 
for c orbitals and, 

a /i \ \ z k z c k x ky 

A mm , (k) = A mm , sin — cos — cos -± 

a f • . ky • k%C 

+ A m m ' sm y sm y sm ~y 



(•5) 



(6) 



A T . k x ky , \ y • ky ky \ k z C 

mm' sm Y C0S 2 mm ' Sm 2 ° 0S 2 / ° 0S T~ 



for m, ml = a, b. Note that beyond the usual p-wave 
symmetry of the sinfe^ and sinfcj, type for the c orbitals, 
we include all three additional p-wave symmetries of the 
sin k/2 type which are induced by the effective attractive 
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interactions between carriers on the neighboring out-of- 
plane Ru orbitals. These interactions are also respon- 
sible for the f-wavc symmetry order parameters, A^ m ,, 
transforming as B\ u in the notation of p^j . This latter 
is symmetry distinct from all p-wave order parameters 
in a tetragonal crystal, unlike some other f-wave states 
which have been proposed ^] . The p z order parame- 
ters A z mml are of A^ u symmetry. In contrast the pairs 
A^ m , , A v mm , for to, to' = a,b are of the same E u sym- 
metry as A^ C ,A^ C . In general the order parameters in 
each distinct irreducible representations will have differ- 
ent transition temperatures. 




b): 











P-band 










p+f 



FIG. 2. Lowest energy eigenvalues, E v (k) on the Fermi sur- 
face; (a) 7 sheet in the plane k z — 0, S 7>ma:c (kp) = 0.25meV, 
i?7,min(kF) = 0.064meV (b) /3 sheet in the plane k x = k y , 
£7/3, max (kir) = 0.32meV. The /—wave order parameter lifts 
the p-wave line nodes (dotted lines) 

Fig. [I] shows the calculated order parameters as a func- 
tion of temperature. One can see that the /—wave and 
p z states have a much smaller T c than the leading p x ,Py 
gap parameters for aa, bb and cc orbitals. Above the 
/-wave T c the order parameters have the symmetries 
A y cc = iA^ c , A bb — iA^ a as expected for a pairing sym- 
metry Q (k x + ik y )e z corresponding to the same time re- 
versal broken pairing state as 3 He — A. The off-diagonal 
components, such as A^ b are small but non-zero, as are 
A bb and A y a . Note that the k-space pairing potentials 
A mm i (k) do not directly correspond to the energy gaps 
on the Fermi surface sheets, shown in Fig. ||, because the 
tight-binding Hamiltonian is non-diagonal in the orbital 
indices. 
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FIG. 3. Calculated specific heat, compared to the experi- 
mental data of NishiZaki et al. .[14] 

Using the above quasi-particle spectra we have calcu- 
lated the specific heat, shown in Fig. ^. Remarkably, al- 
though we fitted to T c only, the specific heat jump at T c 
is 27mJ/mol, in essentially exact agreement with the ex- 
periment |i~4| ]. Thus, contrary to the argument given by 
Zhitomirsky and Rice [[l2| it is not necessary to assume 
a gap opening up on the active 7 sheet only in order 
to obtain the correct discontinuity at T c . Furthermore 
the calculated specific heat also follows the experiment 
closely at lower temperatures, except for the second tran- 
sition at about 0.3T C . Clearly from Fig. [I] this is due to 
the f-wave and p z pairing, which are symmetry decoupled 
from the p x ,p y states, and is not consistent with the ex- 
periments [jl4j. However, if we suppress the /—wave and 
p z pairing then the specific heat agrees perfectly down to 
the lowest temperatures. We believe that this is physi- 
cally reasonable, because the f-wave states will be more 
strongly influenced by impurities than the p-wave states, 
due to a (21 + 1) prefactor in the impurity pair breaking 
parameter 19 1. As can be seen in Fig. |l| the x,y aa and 
cc order paramameters are only slightly affected by the 
presence of the f-wave gap. 

The low temperature limit of the specific heat is power 
law, because our gap parameters have line nodes. These 
are horizontal circles around the cylindrical a and (3 
Fermi surface sheets, as illustrated in Fig. g, while the 
7 sheet is node-less. As can be seen, when the f-wave 
order parameter also becomes finite, the line nodes dis- 
appear. The fact that the slope of specific heat at low 
temperatures (without f) agrees quantitatively with the 
experiments suggests that these horizontal nodes are only 
present on a, /3, as in our model. 

A further independent test of our model is the calcula- 
tion of the superfluid density [Q , shown in Fig. ||. Again 
there is excellent agreement over the whole temperature 
range between T c and zero. Some physical insight into 
the different contributions to the superfluid density can 
be found by setting A^ m , to zero for to = to' = c or 
toto' = a, b and repeating the calculation. The c orbital 
only contribution gives about 70% of the zero tempera- 
ture superfluid density, with the remaining 30% derives 
from the a, b orbital order parameters. One can see that 
the finite slope at zero temperature derives only from a, b 
contributions, consistent with the the line nodes on the 
a, (3 Fermi surfaces. The only disagreement is that the 
absolute magnitude of our calculated zero temperature 
x-y plane penetration depth is only 450A compared to 
1900A determined experimentally pO| , pi]] . This discrep- 
ancy may be due to impurities (the experimental samples 
had T c ~ 1.1 — l.ZK ) or to non-local clcctrodynamic ef- 
fects associated with the line nodes p2[ . We have also 
calculated the temperature dependent zero-field thermal 
conductivity, which is also in good qualitative agreement 
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with the experiments of Izawa et al. ]2 




T/Tc 

FIG. 4. Superfluid density as a function of T (solid line), 
and experimental points from Sample 1 of Bonalde et al. [15] . 
The relative contributions of c and a, b order parameters are 
indicated (dashed lines). 

In summary we would like to emphasize two points. 
Firstly, we have proposed an alternative to the the 'intra- 
band proximity effect' model of Zhitomirsky and Rice JO] 
for describing horizontal line nodes on the a, (3 sheets 
of the Fermi Surface in superconducting Sr 2 Ru04. Our 
bond model differs from theirs in the way the inter- 
layer coupling is implemented. Our description is a real 
space, two point tight-binding interaction such as natu- 
rally arises in any multi-band, extended, negative U Hub- 
bard model, Eq. |l|. To be quite clear about this matter 
we recall that a generic pair-wise interaction like U(r, r'), 
when expressed in the language of a tight-binding model 
Hamiltonian will, in general, give rise to four point inter- 
action parameters Uij^i- The original Hubbard Hamilto- 
nian makes use of the one point parameters Tj\ — Va,n 
whilst the extended Hubbard model is based on two point 

(2) 

parameters E/j ■ — Uij^j. Evidently our 'bond' model is 
a negative U- version of the latter (24| . On the other hand 
the 'proximity effect model' |12j corresponds to using 

(3) — 

Ufj j . The physics of this is often referred to as assisted 
hopping pEfl . If one assumes, as is normally the case in an 
isotropic substance, that U {1) > C/ (2) > C/ (3) > C/ (4) then 
the 'bonds' represent stronger coupling than assisted 
hopping and should be the preferred coupling mecha- 
nism. However, for the tetragonal arrangement of Ru 
atoms in Si~2Ru04 this is no more than a suggestion at 
present. Thus the relative merits of the two models will 
eventually be settled by appeal to experiments. 

Secondly, we wish to stress that in our 'bond' approach 
to the problem the parameters which describe the nor- 
mal state are determined by fitting to the very accurately 
known Fermi Surface and the measured T c determines 
both coupling constants Uu — 40meV, U± = 48meV. 
Thus, the calculated specific heat, Fig. [| and superfluid 
density, Fig. ^, are parameter free quantitative predic- 
tions of the theory. Consequently, their good agreement 
with experiments can be construed as strong support for 



the physical picture represented by our model. Interest- 
ingly a significant feature of this is that the amplitude 
of the gap function, and correspondingly the superfluid 
density, on the a, (3 sheets of the Fermi Surface are com- 
parable to that on the 7 sheet. This is unlike the case 
suggested by the 'inter-band proximity effect' picture (T^j 
where 7 is the active band and a, (3 play a passive role. 
Hopefully, experiments will soon clarify which of the two 
situations prevails. As for the physical mechanism of 
pairing, the fact that U±_ ~ U\\ implies that the pair- 
ing interaction is fairly isotropic in spite of the layered 
structure. 
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